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We consider the one-parameter family of interval maps arising from 
generalized continued fraction expansions known as a-continued fractions. 
For such maps, we perform a numerical study of the behaviour of metric 
entropy as a function of the parameter. The behaviour of entropy is known 
to be quite regular for parameters for which a matching condition on the 
orbits of the endpoints holds. We give a detailed description of the set 
A4 where this condition is met: it consists of a countable union of open 
intervals, corresponding to different combinatorial data, which appear to 
be arranged in a hierarchical structure. Our experimental data suggest 
that the complement of A4 is a proper subset of the set of bounded-type 
numbers, hence it has measure zero. Furthermore, we give evidence that 
the entropy on matching intervals is smooth; on the other hand, we can 
construct points outside of M on which it is not even locally monotone. 

1 Introduction 

Let a£ [0, 1]. We will study the one-parameter family of one-dimensional maps 
of the interval 
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with aj jQ £ N, 6i iQ € {±1} which we call a-continued fraction. These sys- 
tems were introduced by Nakada and are also known in the literature as 
Japanese continued fractions. 



The algorithm, analogously to Gauss' map in the classical case, provides 
rational approximations of real numbers. The convergents ^ =llSL are given by 

f P-l,a = 1 P0,a — Pn+l,a = £n+l,aPn-l,a + a n+l,aPn,a 
\ 9-1, a = Qa.a = 1 Qn+l.a = ^n+l,a1n-l,a + 0-n+l,aQn,a 

It is known (see [9]) that for each a € (0, 1] there exists a unique invariant 
measure fj, a (dx) = p a {x)dx absolutely continuous w.r.t. Lebesgue measure. 




In this paper we will focus on the metric entropy of the T Q 's, which by 
Rohlin's formula ( 13J) is given by 

h{T a ) = -2 log\x\p a (x)dx 

Ja-l 

Equivalcntly, entropy can be thought of as the average exponential growth rate 
of the denominators of convergents: for ;U a -a.e. x € [a — 1, a], 

h(T a ) = 2 lim - log q n . a (x) 

n — *oo TL 

The exact value of h(T a ) has been computed for a > \ by Nakada ([TT]) 
and for \[2 — 1 < a < \ by Cassa, Marmi and Moussa f[TU]1. 

In [9 j , Luzzi and Marmi computed numerically the entropy for a < \[2 — 1 
by approximating the integral in Rohlin's formula with Birkhoff averages 

2 Ar_1 

h(a ) N,x) = --^log\Ti(x)\ 
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for a large number M of starting points x G {a — 1, a) and then averaging over 
the samples: 

1 M 

h(a,N,M) = — ^ft(a,n,a:fc) 
fc=i 

Their computations show a rich structure for the behaviour of the entropy as 
a function of a; it seems that the function a t— > h(T a ) is piecewise regular and 
changes monotonicity on different intervals of regularity. 

These features have been confirmed by some results by Nakada and Natsui 
(\12\, thm. 2) which give a matching condition on the orbits of a and a — 1 

T^(a) = T* 2 (a - 1) for some ki,k 2 G N 

which allows to find countable families of intervals where the entropy is in- 
creasing, decreasing or constant (see section |3J. It is not difficult to check that 
the numerical data computed via Birkhoff theorem fit extremely well with the 
matching intervals of [T2j . 




Figure 2: Numerical data vs. matching intervals 

In this paper we will study the matching condition in great detail. First of all, 
we analyze the mechanism which produces it from a group-theoretical point of 
view and find an algorithm to relate the a-continued fraction expansion of a and 
a — 1 when a matching occurs. This allows us to understand the combinatorics 
behind the matchings once and for all, without having to resort to specific matrix 
identities. As an example, we will explicitly construct a family of matching 
intervals which accumulate on a point different from 0. In fact we also have 
numerical evidence that there exist positive values, such as [0,3,1], which are 
cluster point for intervals of all the three matching types: with k\ < k 2 , k\ = k 2 
and ki > k 2 . 
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We then describe an algorithm to produce a huge quantity of matching 
intervals, whose exact endpoints can be found computationally, and analyze the 
data thus obtained. These data show matching intervals are organized in a 
hierarchical structure, and we will describe a recursive procedure which should 
produce such structure. 

Let now M. be the union of all matching intervals. It has been conjectured 
(C2i sec t- 4, pg. 1213) that M. is an open, dense set of full Lebesgue measure. 
In fact, the correctness of our scheme would imply the following stronger 

Conjecture 1.1. For any n, all elements o/(^rj) have regular continued 

fraction expansion bounded by n. 

Since the set of numbers with bounded continued fraction expansion has 
Lebesgue measure zero, this clearly implies the previous conjecture. 

We will then discuss some consequences of these matchings on the shape 
of the entropy function, coming from a formula in |12j . This formula allows 
us to recover the behaviour of entropy in a neighbourhood of points where a 
matching condition is present. First of all, we will use it to prove that entropy 
has one-sided derivatives at every point belonging to some matching interval, 
and also to recover the exact value of h(T a ) for a > 2/5. In general, though, to 
reconstruct the entropy one also has to know the invariant density at one point. 

As an example, we shall examine the entropy on an interval J on which (by 
previous experiments, see [9], sect. 3) it was thought to be linearly increasing: 
we numerically compute the invariant density for a single value of a € J and 
use it to predict the analytical form of the entropy on J, which in fact happens 
to be not linear. The data produced with this extrapolation method agree with 
high precision, and much better than any linear fit, with the values of h(T a ) 
computed via Birkhoff averages. 

The paper is structured as follows: in section 2 we will discuss numerical 
simulations of the entropy and provide some theoretical framework to justify the 
results; in section 3 we shall analyze the mechanisms which produce the match- 
ing intervals and in section 4 we will numerically produce them and study their 
hierarchical structure; in section 5 we will see how these matching conditions 
affect the entropy function. 
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2 Numerical computation of the entropy 



Let us examine more closely the algorithm used in [9 to compute the entropy. 
A numerical problem in evaluating Birkhoff averages arises from the fact that 
the orbit of a point can fall very close to the origin: the computer will not 
distinguish a very small value from zero. In this case we neglect this point, and 
complete the (pseudo)orbit restarting from a new random seecj^J As a matter 
of fact this algorithm produces an approximate value of 



h e (a) := / f e (x)d^ a (x) with f e (x) 



|x| < € 

-2 log |x| jxj > e 

where e = 1CP 16 ; of course h e (a) is an excellent approximation of the entropy 
h(a), since the difference is of order eloge -1 . To calculate h e (a) we use the 
Birkhoff sums 

1 N ~ 1 

h e (a,N,x) := - £ MTi(x)) 

3=0 

and in [16] the fourth author proves that for large N the random variable 
h(e,N,-) is distributed around its mean h e (a) approximately with normal law 
and standard deviation a e (a) / yN where 

lim / f ^'-"/^ ) 2 ^ 

n->+ooJ Ia \ y/n J 

which explains the aforementioned result by Luzzi and Marmi 9 . 

One of our goals is to study the function a i— > of (a), in particular we ask 
whether it displays some regularity like continuity or semicontinuity. To this 
aim we pushed the same scheme as in [9] to get higher precision: 

1. We take a sample of values a chosen in a particular subinterval J C [0, 1]; 

2. For each value a we choose a random sample {x\, ...,xm} in la (the car- 
dinality M of this sample is usually 10 6 or 10 7 ) ; 

3. For each n € I a (i= 1, M) we evaluate h e (a, N, x t ) as described before 
(the number of iterates N will be 10 4 ); 

4. Finally, we evaluate the (approximate) entropy and take record of standard 
deviation as well: 



1 M 

h e (a, N, M) :— — ^ h e (a, N, x t ) 



i=i 



a t (a) 



\ 



1 M 



i=i 



2 Anothor choice is to throw away the whole orbit and restart; it seems there is not much 
difference on the final result 



5 



2.1 Central limit theorem 



Let us restate more precisely the convergence result for Birkhoff sums proved in 
[IB] . Let us denote by BV(I a ) the space of real-valued, /i Q -integrable, bounded 
variation functions of the interval I a . We will denote by S n f the Birkhoff sum 

n-l 
3=0 

Lemma 2.1. Let a S (0, 1] and f be an element of BV(I a ). Then the sequence 

M n = l ( M^V^ )V 

converges to a real nonnegative value, which will be denoted by a 2 . Moreover, 
a 2 = if and only if there exists u G L 2 (fi a ) such that up a £ BV{I a ) and 



f- 



fd-Ho 



U — U O T n 



(1) 



The condition given by (JT|) is the same as in the proof of the central limit 
theorem for Gauss' map, and it's known as cohomological equation. The main 
statement of the theorem is the following: 

Theorem 2.2. Let a 6 (0, 1] and f be an element of BV(L a ) such that has 
no solutions. Then, for every v € K we have 



lim /j, a 



Snf -nfffd/ia 



< V 



'2-k J_ 



- x *' 2 dx 



Since we know that the invariant density p a is bounded from below by a 
nonzero constant, we can show that 

Proposition 2.3. For every real-valued nonconstant f € BV(L a ), the equation 
has no solutions. Hence, the central limit theorem holds. 

Now, for every e > the function f e define in the previous section is of 
bounded variation, hence the central limit theorem holds and the distribution 
of the approximate entropy h t (a,N,-) approaches a Gaussian when N — > oo. 
As a corollary, for the standard deviation of Birkhoff averages 



Std 


Snfe 


= E 


( n Sj^ a ) 




n 







1 1/2 



1 



2.2 Speed of convergence 

In terms of numerical simulations it is of primary importance to estimate the 
difference between the sum computed at the n th step and the asymptotic value: 
a semi-explicit bound is given by the following 
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Theorem 2.4. For every nonconstant real-valued f G BV(I a ), there exists 
C > smc/i that 



sup 



S n f -n J j fd[j, c 



< V 



1 



2 dx 



< 



C 



Proof. It follows from a Berry-Esseen type of inequality. For details see ([T], 
th.8.1). □ 



2.3 Dependence of standard deviation on a 

Given these convergence results for the entropy, it is natural to ask how the 
standard deviation varies with a. In this case not a single exact value of <7 £ (a) 
is known; by using the fact that natural extensions of T a are conjugate ([7], 
[12]). it is straightforward to prove the 

Lemma 2.5. The map a t— * o~(a) is constant for a G [V2— 1, v/ ^~ 1 ] ■ 

Proof. See appendix. □ 

The numerical study of this quantity is pretty interesting. We first consid- 
ered the window J — [0.295, 0.304299], where the entropy is non-monotone. On 
this interval the standard deviation shows quite a strange behaviour: the values 
we have recorded do not form a cloud clustering around a continuous line (like 
for the entropy) but they cluster all above it. 

One might guess that this is due to the fact that the map a i— > a (a) is 
only semicontinuous, but the same kind of asymmetry takes place also on the 
interval J = [0.616,0.618], where a 2 is constant. Indeed, we can observe the 
same phenomenon also evaluating <7 e (a) for a fixed value a but taking several 
different sample sets. 

On the other hand this strange behaviour cannot be detected for other maps, 
like the logistic map, and could yet not be explained. Nevertheless, we point 
out that if you only consider C 1 observables, the standard deviation of Birkhoff 
sums can be proved continuous, at least for a € (0.056,2/3). See [IB"] . 



3 Matching conditions 

In [T5], Nakada and Natsui found a condition on the orbits of a and a — 1 
which allows one to predict more precisely the behaviour of the entropy. Let 
us denote for any a G [0, 1], x G I a , n > 1 by M aiXtn the matrix such that 
T£(x) = M-X n (x), i.e. 

They proved the following 
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Figure 3: Variance on the interval J = [0.295,0.304299]. 

Theorem 3.1. ( [12j. thm. 2) Let us suppose that there exist positive integers 
k\ and k 2 such that 

(I) {T£(a) : < n < ki} n {T™(a - 1) : < m < k 2 } = 

(II) M a<aM = ( J } ) M a , a _ life2 [=► T**(a) = T^(a - 1) ] 

(III) T^(a) [= T*»(a-1)] £ {a,a-l} 

T/ien f/iere exists rj > smc/i £/iai, on (a — 77, a + 77), h(T a ) is : 

(i) strictly increasing if k± < k 2 

(ii) constant if k\ = k 2 

(iii) strictly decreasing if k\ > k 2 

It turns out that conditions (I)- (II)- (III) define a collection of open intervals 
(called matching intervals); they also proved that each of the cases (i), (ii) and 
(iii) takes place at least on one infinite family of disjoint matching intervals clus- 
tering at the origin, thus proving the non-monotonicity of the entropy function. 
Moreover, they conjectured that the union of all matching intervals is a dense, 
open subset of [0, 1] with full Lebesgue measure. 
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Figure 4: Plot of the standard deviation of the different runs on the Gauss map 



In the following wc will analyze more closely the mechanism which leads to 
the existence of such matchings. As a consequence, we shall see that it looks 
more natural to drop condition (III) from the previous definition and replace 
(II) with 

(II') M aiClM -l=±( 1 Q \ ) Ma.a-x,*,-! 

(which implies Tkl li (a) + T^k^i) = 

We can now define the matching set as 

M = {aE (0, 1] s. t. (I) and (II') hold } 

Note M. is open, since the symbolic codings of a up to step k\ — 1 and of a — 1 
up to step ki — 1 are locally constant. 

Moreover, we will see that under this condition it is possible to predict the 
symbolic orbit of a — 1 given the symbolic orbit of a, and viceversa. As an 
application, we will construct a countable family of intervals which accumulates 
in a point different from 0. 

Let us point out that our definition of matching produces a set slightly bigger 
than the union of all matching intervals satisfying condition (I, II, III): in fact the 
difference is just a countable set of points. 
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3.1 Structure of PGL(2, Z) 

Let us define PGL(2,Z) := GL(2, 1)/{±I}, PSL(2,Z) := SL(2,Z)/{±I}. We 
have an exact sequence 

1 -> PSL(2,Z) -» PGL(2,Z) -> {±1} -> 1 

where the first arrow is the inclusion and the second the determinant; moreover, 
if we consider the group 

c-{( ^ MU z x z 

° ~ H ±1 J / - 2Z X 2Z 
and let C = C/{±7}, then 

Cn?SL(2,Z) = {e} 
therefore we have the semidirect product decomposition 

PGL(2,Z) = PSL(2,Z) x C 
Now, it is well known that PSL(2,Z) is the free product 

PSL(2,Z) =< S > * <U > 

where 




are such that S 2 = I, U 3 = I. Geometrically, S represents the function {z — > 
— -}, and if we denote by T the element corresponding to the translation {z — > 
z+1}, we have U = ST. 

The matrix V = ^ ^ ^ ^ projects to a generator of C and it satisfies 

V 2 = I, VSV- 1 = VSV = S and VTV -1 = T^ 1 in PGL(2,Z) so we get the 
presentation 

PGL(2,Z) = {S, T, V | S 2 = 7, (ST) 3 = /, V 2 = /,VSV _1 = S, VTV -1 = T" 1 } 



3.2 Encoding of matchings 

Every step of the algorithm generating a-continued fractions consists of an op- 
eration of the type: 

z - - c e e {±l},c e N 

z 

which corresponds to the matrix T~ c SV e ^ with 

. , f if e = -1 
6 ^= 1 if e = l 
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so if x belongs to the cylinder ((ci, e\), . . . , (cfc, e&)) we can express 

T*(x) = T- Ck SV e{tk) . ..T~ cl SV e( - ei) (x) 

Now, suppose we have a matching T^ 1 (a) — T^ 2 (a — 1) and let a be- 
long to the cylinder ((ai,ei), . . . , (a^e/ci)) ano - ol — 1 belong to the cylinder 
((61, . . . , (bk 2 ,T]k 2 )). One can rewrite the matching condition as 

rp-a hl gye(e hl ) _ _ _ j.-ai gye^i) ( a \ _ T~bk 2 gye(vk 2 ) _ _ _ ji-Jj gye( m )y-l 

hence it is sufficient to have an equality of the two Mobius transformations 

rp-a-kx gy^kx) _ _ _ rp-a-i gye(ei) _ y-bfc 2 gye(nk 2 ) _ _ _ 5^(771) 

We call such a matching an algebraic matching. Now, numerical evidence shows 
that, if a matching occurs, then 

ei = +1 

e* = -1 for 2 < i < fci - 1 
r/i = — 1 for 1 < i < k2 — 1 

If we make this assumption we can rewrite the matching condition as 

_ y e (v k2 )rpb k2 (_i)' e(l,fc 2 , + gy-6fc 2 -i o . . . ST^ 1 

which implies e(efc 1 ) = e(j]k 2 ) + 1, i.e. £k 1 ilk 2 = ~ 1- If f° r instance e(efe 1 ) = 1 
and e(?7fe 2 ) = 0, by substituting T = SC/ one has 

(u 2 s) a »i u{su) ak ^- 2 su 2 . . . su 2 (su) ai - 2 sus = 
= {u 2 s) bh *- l us{u 2 s) bk *- i - 2 us . . . US(U 2 S) bl - 2 US 

Since every element of PSL{2, Z) can be written as a product of S and Z7 in 
a unique way, one can get a relation between the a r and & r . Notice that, since 
we are interested in a < \ 2— 1, Oj > 2 and bi > 2 for every i, hence there is no 
cancellation in the equation above. By counting the number of (U 2 S) blocks at 
the beginning of the word, one has = b^ 2 — 1, and by semplifying, 

{suy*i-i- 2 su 2 . . . su 2 (su) ai - 2 sus = , . 

= 5(C/ 2 5) bfe 2-i- 2 C/5 . . . US{U 2 S)^- 2 US { > 

If one has e(ej ;i ) = and e(?7fc 2 ) = 1 instead, the matching condition is 

(SU) ak i- l SU 2 {SU) ak ^- 2 SU 2 . . . SU 2 (SU) ai - 2 SUS = 

(SU) bk ? stfs^sf^-^us . . . US(U 2 S) bl - 2 US 

which implies bk 2 — &&1 — 1, and simplifying still yields equation ([2]). 
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Let us remark that ([2]) is equivalent to 

T -l ST -a H - lS _ _ _ y-ai^y = VST^^S. . . T~ bl ST' 1 

which is precisely condition (II'): by evaluating both sides on a 

1 1 

r*i-i(a) + T^-^a- 1) ~ 

Moreover, from |2]) one has that to every a r bigger than 2 it corresponds 
exactly a sequence of bi — 2 of length precisely a r — 2, and viceversa. More 
formally, one can give the following algorithm to produce the coding of the 
orbit of a — 1 up to step — 1 given the coding of the orbit of a up to step 
k\ — 1 (under the hypothesis that an algebraic matching occurs, and at least k\ 
is known). 

1. Write down the coding of a from step 1 to k\ — 1, separated by a symbol * 



2. Subtract 2 from every a r ; if a r = 2, then leave the space empty instead of 
writing 0. 

ai — 2*a2 — 2 * ■ ■ ■ * akx-i — 2 

3. Replace stars with numbers and viceversa (replace the number n with n 
consecutive stars, and write the number n in place of n stars in a row) 

4. Add 2 to every number you find and remove the stars: you'll get the 
sequence . . . , fo^ 3 _i) . 

Example. Let us suppose there is a matching with k\ = n + 3 and a has initial 
coding ((3, +), (4, — )", (2, — )). The steps of the algorithm are: 

Step 1 

3*4*4*---*4*2 



Step 2 



Step 3 

★1** 1 **1 . . . 1 **1 

V v ' 

n times 

Step 4 

2 3 2 3 ...2 3 



so 



the coding of a — 1 is ((2, — )(3, — )) n+1 , and k^ = 2n + 3. 
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3.3 Construction of matchings 

Let us now use this knowledge to construct explicitly an infinite family of match- 
ing intervals which accumulates on a non-zero value of a. For every n, let us con- 
sider the values of a such that a belongs to the cylinder ((3, +), (4, — )™, (2, — )) 
with the respect to T a . Let us compute the endpoints of such a cylinder. 



The right endpoint is defined by 

-4 -1 
1 



-3 1 
1 



i.e. 



1 1 
1 



-4 -1 
1 



(a) = a — 1 

(a) = a 



The left endpoint is defined by 

-4 -1 
1 



-3 1 
1 



i.e. 



-2 -1 
1 



-4 -1 
1 



-3 1 
1 



(«) = - 



-3 1 
1 



a + 2 



(a) = a 



By diagonalizing the matrices and computing the powers one can compute 
these value explicitly. In particular, 



a, 



73 - 1 40\/3 - 69 

2 + 13 
\/3-l 10\/3-12 



(2 + V3)- 2n + 0((2 + V3)- 4n ) 
(2 + V3)- 2n + 0((2 + V3)- in ) 



2 13 

The as such that a— I belongs to the cylinder ((2, — ), (3, — ))" +1 are defined 
by the equations 



-3 -1 
1 

for the left endpoint and 



-3 -1 
1 



-2 -1 
1 



-2 -1 
1 



-i 71+1 



(a- 1) = a- 1 



71+1 



(a — 1) = a 



for the right endpoint, so the left endpoint corresponds to the periodic point 
such that 



-3 -1 
1 



-2 -1 
1 



(a- 1) = a - 1 



i.e. 



V3-1 



13 



and 



o? max = ^ + ^i^(2 + V3)-^ + 0((2 + ^r 4 ") 
By comparing the first order terms one gets asymptotically 

a 2 < a 1 < o? < a 1 

hence the two intervals intersect for infinitely many n, producing infinitely many 
matching intervals which accumulate at the point a = ^ 2 ~ 1 ■ The length of 
such intervals is 

o? max - «Ln = 567 ' 2 3 6 2?V3 (2 + ^)- 2 " + 0((2 + V3)^) 

4 Numerical production of matchings 

In this section we will describe an algorithm to produce a lot of matching in- 
tervals (i.e. find out their endpoints exactly), as well as the results we obtained 
through its implementation. Our first attempt to find matching intervals used 
the following scheme: 

1. We generate a random seed of values on belonging to [0,1] (or some other 
interval of interest). When a high precision is needed (we manage to 
detect intervals of size 1CP 60 ) the random seed is composed by algebraic 
numbers, in order to allow symbolic (i.e. non floating-point) computation. 

2. We find numerically candidates for the values of k\ and k 2 (if any) simply 
by computing the orbits of a and of a — 1 up to some finite number of steps, 
and numerically checking if T^ 1 (a) = T^ 2 (a — 1) holds approximately for 
some ki and k 2 smaller than some bound. 

3. Given any triplet (a, ki,k 2 ) determined as above, we compute the symbolic 
orbit of a up to step k\ — 1 and the orbit of a — 1 up to step k 2 — 1 . 

4. We check that the two Mobius transformations associated to these sym- 
bolic orbits satisfy condition (II'): 

1 (I 



= ± I M aia _i i / S2 _i 



We solve the system of quadratic equations which correspond to imposing 
that a and a — 1 have the same symbolic orbit as a and a — 1 , respectively. 

Let us remark that this is the heaviest step of the whole procedure since 
we must solve fci + k 2 — 2 quadratic inequalities; for this reason the value 
k = k\ + k 2 may be thought of as a measure of the computational cost of 
the matching interval and will be referred to as order of matching. 
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Following this scheme, we detected more than 10 7 matching intervals, whose 
endpoints are quadratic surds; their union still leaves many gaps, each of which 
smaller than 6.6 • 10~ 6 . A table with a sample of such data is contained in the 
appendix. [^] 

In order to detect some patterns in the data, let us plot the size of these 
intervals (figure [5j. For each matching interval ]a_,a+[, we drew the point of 
coordinates (o!_, a+ — a_). 




0.05 0.06 0.07 0.08 



Figure 5: Size of matchings 



It seems there is some self-similar pattern: in order to understand better 
its structure it is useful to identify some "borderline" families of points. The 
most evident family is the one that appears as the higher line of points in 
the above figure (which we have highlighted in green): these points correspond 
to matching intervals which contain the values 1/n, and their endpoints are 
n], a + (n) ■■ 



a_(n) = \[\fn 2 + 4 — n], a + (n) = _ [y/n 2 + 2n — 3 — n + 1]; this is the 
family I n already exhibited in |12j . Since a_(n) = 1/n — 1/n 3 + o(l/n 3 ) and 
a+(n) = 1/n + 1/n 3 + o(l/n 3 ), for n ^ 1 the points (a_(n), ot+(n) — a_(n)) 
are very close to This suggests that this family will "straighten" if we 

replot our data in log-log scale. This is indeed the case, and in fact it seems 
that there are also other families which get perfectly aligned along parallel lines 
of slope 3 (see figure [6]) . 

If we consider the ordinary continued fraction expansion of the elements of 
these families we realize that they obey to some very simple empiricaQ rules: 

(i) the endpoints of any matching interval have a purely periodic continued 
fraction expansion of the type [0, oi, a,2, a m , 1] and [0, ai, 02, a m + 1]; 



3 A more efficient algorithm, which avoids random sampling, will be discussed in subsec- 
tion [41] 

4 U niortunately we are still not able to prove all these rules. 
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Figure 6: Same picture, in log-log scale. 



this implies that the rational number corresponding to [0, di, ci2, a m + l] 
is a common convergent of both endpoints and is the rational with smallest 
denominator which falls inside the matching interval; 

(ii) any endpoint [0, ai, 02, a m ] of a matching interval belongs to a family 
{[0, a, a2, a m ] : a > max2<i< m ai}; in particular this family has a 
member in each cylinder B n :— {a : l/(n + 1) < a < 1/n} for n > a, so 
that each family will cluster at the origin. 

(ii ! ) other families can be detected in terms of the continued fraction expansion: 
for instance on each cylinder B n (n > 3) the largest matching interval on 
which h is decreasing has endpoints with expansion [0, n, 2, 1, n — 1,1] and 
[0,n,2,l,n] 

(iii) matching intervals seem to be organized in a binary tree structure, which 
is related to the Stern-Brocot tre^] one can thus design a bisection al- 
gorithm to fill in the gaps between intervals, and what it's left over is a 
closed, nowhere dense set. This and the following points will be analyzed 



extensively in subsection 4.1 



(iv) if a € B n is the endpoint of some matching interval then a — [0; a%, 0,2, a„ 
with a, < n Vi £ {1, m}; this would imply that the values a £ B n which 
do not belong to any matching interval must be bounded-type numbers 
with partial quotients bounded above by n; 

(v) it is possible to compute the exponent (fci,/^) of a matching from the 
continued fraction expansion of any one of its endpoints. 



5 Sometimes also known as Farey tree. See [3]. 
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Figure 7: A few of the largest matching intervals in the window [1/4, 1/3], and the 
corresponding nodes of Stern-Brocot tree. The value on the y axis is the logarithm of 
the size of each interval. 



From our data it is also evident that the size of these intervals decreases as 
ki + hi increases, and low order matchings tend to disappear as a approaches 
zero. Moreover, as a tends to the space covered by members of "old" families 
of type (ii) encountered decreases, hence new families have to appear. One can 
quantify this phenomenon from figure [6] since the size of matching intervals in 
any family decreases as 1/n 3 on the interval cylinder B n (whose size decreases 
like 1/n 2 ): this means that, as n increases, the mass of B n gets more and more 
split among a huge number of tiny intervals. 

This fact compromises our numerical algorithm: it is clear that choosing 
floating point values at random becomes a hopeless strategy when approaching 
zero. Indeed, even if there still are intervals bigger than the double-precision 
threshold, in most cases the random seed will fall in a really tiny interval corre- 
sponding to a very high matching order: this amounts to having very little gain 
as the result of a really heavy computation. 

We still can try to test numerically the conjecture that the matching set has 
full measure on [0, 1]; but we must expect that the percentage of space covered 
by matching intervals (found numerically) will decrease dramatically near the 
origin, since we only detect intervals with k\ + ki bounded by some threshold. 
The matching intervals we have found so far cover a portion of 0.884 of the 
interval [0, 1]; this ratio increases to 0.989 if we restrict to the interval [0.1, 1] 
and it reaches 0.9989 restricting to the interval [0.2, 1]. 
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Figure 8: Dependence of the order k = k\ + k% of a matching interval on the left 
endpoint 

The following graph represents the percentage of the interval [x, 1] which is 
covered by matching intervals of order k = k\ + &2 for different values of It 
gives an idea of the gain, in terms of the total size covered by matching intervals, 
one gets when refining the gaps (i.e. considering matching intervals of higher 
order) . 




Figure 9: Percentage of covering by matching intervals 

Finally, to have a more precise idea of the relationship between order of 
matching and size of the matching interval it is useful to see the following 
scattered plot: the red dots correspond to matching intervals found using a 

6 Let us point out that for big values of k the graph does not take into account all matching 
intervals of order k but only those we have found so far. 
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random seed, and the green ones to intervals found using the bisection algorithm. 
The two lines bounding the cloud correspond to matching intervals with very 
definite patterns: the upper line corresponds to the family /„ (with endpoints 
of type [0; n] and [0; n — 1, 1]), the lower line corresponds to matching intervals 
with endpoints of type [0; 2, 1, 1, 1, 1, 1] and [0; 2, 1, 1, 1, 2]. The latter ones 
converge to 3 ~ 2 , which is the supremum of all values where the entropy is 
increasing. 

Thus numerical evidence shows that, if J is an interval with matching order 
k = fci + &2 then the size of J is bounded below by |J| > coe~ Clk where 
c = 8.4423... and c\ — 0.9624.... On the other hand we know for sure that, 
on the right of 0.0475 (which corresponds to the leftmost matching interval of 
our list), the biggest gap left by the matching intervals found so far is of order 
6.6 • 10~ 6 . So, if J is a matching interval which still does not belong to our list, 
either J C [0, 0.0475] and k > 20 (see figure j8|, or its size must be smaller than 
6.6 • 10~ 6 and by the forementioned empirical rule, its order must be k > 14.6. 
Hence, our list should include all matching intervals with ki + ki < 14. 















^^^^ 
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Figure 10: The order k\ + ki versus the logarithm of the size of the first 10 7 matching 
intervals found. 



4.1 The matching tree 

As mentioned before, it seems that matching intervals are organized in a bi- 
nary tree structure. To describe such structure, we will provide an algorithm 
which allows to construct all matching intervals by recursively "filling the gaps" 
between matching intervals previously obtained, similarly to the way the usual 
Cantor middle third set is constructed. 

In order to do so, let us first notice that every rational value r € Q has two 
(standard) continued fraction expansions: 

r= [0;a 1 ,a 2 ,...,a m ,l] = [0; a x , a 2 , a m + 1] 
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One can associate to r the interval whose endpoints are the two quadratic surds 
with continued fraction obtained by endless repetition of the two expansions of 
r: 

Definition 4.1. Given r € Q with continued fraction expansion as above, we 
define I r to be the interval with endpoints 

[0;ai,a 2 ,...,a TO ,l] and [0; 01,02,..., o m + 1] 

(in any order). The strings Si := {ai, . . . , a m , 1} and S2 ■= {ai, • • • , a m + 1} 
will be said to be conjugate and we will write S2 = (Si)' . 

Notice that r e I r . It looks like all matching intervals are of type I r for some 
rational r. On the other hand, 

Definition 4.2. Given an open interval I 3 [0, 1] one can define the pseudocen- 
ter of I as the rational number r e / n Q which has the minimum denominator 
among all rational numbers contained in I. 

It is straightforward to prove that the pseudocenter of an interval is unique, and 
the pseudocenter of I r is r itself. 

We are now ready to describe the algorithm: 

1. The rightmost matching interval is [ v/ ^~ 1 , 1]; its complement is the gap 

J =[o,4pi]. 

2. Suppose we are given a finite set of intervals, called gaps of level n, so 
that their complement is a union of matching intervals. Given each gap 
J = [a~,a + ], we determine its pseudocenter r. Let cr 1 — [0; S, a ± , S^] 
be the continued fraction expansion of a^, where S is the finite string 
containing the first common partial quotients, a + ^ a~ the first partial 
quotient on which the two values differ, and the rest of the expansion of 
a* , respectively. The pseudocenter of [a~ , a + ] will be the rational number 
r with expansions [0; 5*, a, 1] = p/q = [0; S*, a + 1] where a := min(a + , a~). 

3. We remove from the gap J the matching interval I r corresponding to the 
pseudocenter r: in this way the complement of I r in J will consist of two 
intervals Ji and J2 , which we will add to the list of gaps of level n + 1 . It 
might occur that one of these new intervals consists of only one point, i.e. 
two matching intervals are adjacent. 

By iterating this procedure, after n steps we will get a finite set Q n of gaps, 
and clearly \J Je g n 1 J C \J Je g n J- We conjecture all intervals obtained by 
taking pscudocenters of gaps are matching intervals, and that the set on which 
matching fails is the intersection 

Go° ■= n u j < 



20 



The next table contains the list of the elements of the family Q n of gaps of 
level n for n — 0..4: when a gap is reduced to a point we mark the corresponding 
line with the symbol *. 
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We still cannot prove that this is the right scheme, but the numerical ev- 
idence is quite robust: all 1.1 • 10 6 intervals obtained by running the first 23 
steps, for instance, turn out to be real matching interval^ 

We can also prove the 

Lemma 4.1. Goo consists of numbers of bounded type; more precisely, the ele- 
ments of Goo H (^xj, 4] have regular continued fraction bounded by n. 

Proof. The scheme described before forces all endpoints of matching intervals 
containtcd in the cylinder B n =]l/(rt + l), l/n[ to have quotients bounded by n. 
We now claim that, if 7 = [0; ci, C2, c„, ...] ^ M, then, cj, < ci for all fceN. 

If 7 ^ Ai then 7 € Ujee ^ ^ or a ^ n E W, let us call J„ the member 
of the family Q n containing 7. It may happen that there exists uq such that 
J n = {7} Vrt > no, so. 7 is an endpoint of two adiacent matching intervals, 

7 We compared them to the list obtained as in section Wl 



21 




Figure 11: Recursive construction of the matching set 



hence it has bounded type. Otherwise, J„ = [a n , /3 n ] with f3 n — a n >0 Vn > c\, 
where a n ,f3 n are the endpoints of two matching intervals. Now, \i p n /q n is the 
pseudocenter of J„ from the minimality of q n it follows that \[3 n — a n \ < 2/q n , 
but also that q n +i > In (since Pn+i/ln+i £ Jn+i C J n ); these two properties 
together imply that < 7 — a n < 2/q n — > as n — > +00. This implies 7 
cannot be rational, since 7 S J„ Vn and the minimum denominator of a rational 
sitting in J n is g n — > +00. Hence, since a n — > 7, for every fixed fc G N, there is 
some n(fc) such that for all n > n(k) all the partial quotients up to level k of 7 
coincide with those of a n , which are bounded by c\ . □ 

As a consequenc e, th e validity of our algorithm (Qoo — [0, 1] \ JA) would 
imply the conjecture |l.l| |^] 

Notice Qoo n (l/(n+ 1), 1/n] has Hausdorff dimension strictly smaller than 
one for each n. Moreover, the Hausdorff dimension of n-bounded numbers tends 
to 1 as n — > 00. We think that, similarly, H.dim{(^q-j-, ^] \ A4} — > 1: this would 
explain why finding matching intervals near the origin becomes a tough task. 

Remark. Since we have associated a rational number to each matching interval, 
one can think of the bisection algorithm as acting on Q, and get a binary tree 

8 Our conjecture implies that also the set where the original conditions given by Nakada- 
Natsui hold has full measure; the equivalent of lemma |4.1| is, however, not true for their 
matching set, which differs from ours for a countable number of points. 
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whose nodes are rationals: this object is related to the well-known Stern-Brocot 
tree. (For an introduction to it, see J^). 

Given that all matching intervals correspond to some rational number, one 
can ask which subset of Q actually arises in that way. 

Definition 4.3. An interval I r , r € Q is maximal if I r ~D I r / Vr' € I r n Q. 

We conjecture that the matching intervals are precisely the maximal inter- 
vals, so that the matching set is 

M= (J I r = |J Ir 

re[0,l]nQ re[0,i]nQ 
I r maximal 

As a matter of fact we can actually prove that the complement of the family 
Q n produced by the bisection algorithm consists of a family of maximal intervals: 
the proof of this fact is rather technical and will appear in a forthcoming paper. 

We have also found an empirical rule to reconstruct the periods {k\, kz) of a 
matching interval from the labels of its enpoints. Let S — [ai,...,at] be a label 
of the endpoint s of some matching interval: 

1. If s is a left endpoint then 

fci = 2 + >J aj , fc 2 = dj . 

j even j odd 



2. If s is a right endpoint then 



ki = 1 + Oj, fc 2 = 1 + ( 

j even j odd 



Trusting this rule, we are able to prove that every neighbourhood of the 
point [0,3, 1] contains intervals of matching of all types: with k\ < k^, fci = ki 
and k\ > k%. Indeed, it is not difficult to realize that [0, 3, 1] is contained in the 
family of gaps Jp of endpoints [0, 3, P] and [0, 3, P, 1] where P is a string of the 
type 1, 1, 1, 1 of even length; by our rule the left endpoint of Jp is the right 
endpoint of an interval of matching where k\ < ki. Nevertheless, performing a 
few steps of the algorithm, it is not difficult to check that the gap Jp contains 
the interval Cp of enpoints [0, 3, P, 2, 1, 1] and [0, 3, P, 2, 1, 1] (on which k% = hz) 
but also Dp of enpoints [0, 3, P, 2, 1, 2, l]and [0, 3, P, 2, 1, 3] (on which k\ > £2). 



4.2 Adjacent intervals and period doubling 

Let us now focus on pairs of adjacent intervals (corresponding to isolated points 
in [0, 1]\.M): our data show they all come in infinite chains, and can be obtained 
from some starting matching interval via a "period doubling" construction. 
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Let's start with a matching interval ]a, f3[ ; a = [0; S) where S is a sequence of 
positive integer of odd length; define the sequence of strings 

\ Sq = S 

\ <Sn+l = {S n S n ) 

where S' denotes the conjugate of S as in def. 
b n := [0; S' n ]; then the sequence /„ :=]a n ,b n [ is formed by a chain of adjacent 
intervals: clearly b n +i = a n , moreover a n < b n because \S n \ is odd for all n. 

Assuming this scheme, we can construct many cluster points of matching 
intervals. For instance, let us look at the first (i.e. rightmost) one: we start 
with the interval ](v5 — l)/2, 1[ so that the first terms of the sequence S n are 

50 = (1) 

51 = (2) 

5 3 = (2,1,1) 

5 4 = (2,1,1,2,2) 

5 5 = (2,1,1,2,2,2,1,1,2,1,1) 

The corresponding sequence a n converges to the first (i.e. rightmost) point 
a where intervals of matching cluster. We can also determine the continued 
fraction expansion of the value a, since it can be obtained just merging the 
strings (S n ) n S N 

a = [0,2,1,1,2,2,2,1,1,2,1,1,2,1,1,2,2,2,1,1,2,2,2,1,1,2,2,2,1,1,2,1,1,2,1,1,2,2,2,...] 

Numerically a ^ 0.386749970714300706171524803485580939661. . . 

It is evident from formula ^ that any such cluster point will be a bounded- 
type number; one can indeed prove that no cluster point of this type is a 
quadratic surd. 



4.1 



Let a n :— [0; S n ] and 



5 Behaviour of entropy inside the matching set 

In [12., the following formula is used to relate the change of entropy between 
two sufficiently close values of a to the invariant measure corresponding to one 
of these values: more precisely 



Proposition 5.1. Let us suppose the hypotheses of prop. \3.1\ hold for a: then 
for r\ > small enough 



h(T a ) 

1 + (k 2 - ki)n a ([a ~ i], a}) 
and similarly 

h(T ) = h(T a+ri ) 

1 al l + te-fcO/Wk « + »?]) 

By exploiting these formulas, we will get some results on the behaviour of 
h(T a ). 

9 This can be done since, by S n is a substring of S n +i- 

10 This pattern has been checked up to level 10, which corres pon ds to a matching interval 
of size smaller than 10 -200 ; see also the second table in section 6.1 
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5.1 One-sided differentiability of h(T a ) 

Equation Q has interesting consequences on the differentability of h: we can 
rewrite it as 

h(T a ) - tyTa-r,) = h(T a - v )(k 2 - h)p, a ([a - r?, a}) 
and dividing by i] 

h(T a ) h(T a -„) = ^([a-nM) 

Tj Tj 

Since p a has bounded variation, then there exists R(a) — lim x ^ a - p a (x), there- 
fore 

lim = R{a) 

and by the continuity of h (which is obvious in this case by equation ([4] 
lim h (T a )-h(T a ^) =h{Ta){k2 _ ki) Um pa[x) 

hence the function a i— > h(T a ) is left differcntiable in a. On the other hand, one 
can slightly modify the proof of ^ and realize it is equivalent to 

h(T a+v ) — 



1 + (fci - k 2 )p a {[a - 1, a - 1 + rj\) 
which reduces to 

h(T a+v ) - h(T a ) _ fi a ([a - l,a - 1 + rj\) h(T a )(k 2 - fci) 

Tj rj 1 + (fci - k 2 )p a ([a - l,a - 1 + rj\) 

Since the limit 

fi a ([a - 1, a - 1 + r/}) 

lim = lim Pala;) 

»7->0 r/ •£_>(<*_!)+ 

also exists, then h(T a ) is also right differcntiable in a, more precisely 
lim /t(T °" H ' ) ~ KTa) = h(T a )(k 2 - h) lim p a (x) 

77— >0 7^ a;— + 

We conjecture that in such points the left and right derivatives are equal. 
This is trivial for k\ = k 2 ; for k\ 7^ k 2 it is equivalent to say lim c _ K! ,- p a [x) = 

]im x -,( a -l)+ Pa(x). 
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5.2 The entropy for a > | 

Corollary 5.2. For | < a < \/2 — 1, £/ie entropy is 

fc(T a ) = 



61og(4±!) 



Proof. Every a in the interval (0.4, \/2 — 1) satisfies the hypotheses of the the- 
orem with fci = fc 2 = 3, hence h(T a ) is locally constant, and by continuity 
h(T a ) = h(T^_ 1 ), whose value was already known. 

□ 

Remark. By using our computer-generated matching intervals, we can analo- 
gously prove h(T a ) = h{T^_ x ) for^-l>a> 0.386749970714300706171524... 

5.3 Invariant densities 

In the case a > \/2 — 1 it is known that invariant densities are of the form 

r 
i— 1 

where the Ii are subintervals of [a — 1 , a] . 

For these values of a, a matching condition is present and the endpoints of 
the Ii (i.e. the values where the density may "jump") correspond exactly to 
the first few iterates of a and a — 1 under the action of T a . We present some 
numerical evidence in order to support the 

Conjecture 5.3. Let a G [0, 1] be a value such that one has a matching of type 
(fci,fe) (i-e. with T^ 1 (a) — T* 2 (a— \)). Then the invariant density has the 
form 

Pa{x) = ^Xh{x) (6) 

i— 1 

where each Ii is an interval with endpoints contained in the set 

S := {T™(a) : < m < fci} U {T£(a - 1) : < n < k 2 } 
Therefore, the number of branches is bounded above by k\ + k 2 — 1 . 

In all known cases, moreover, there exists exactly one Ii which contains a 
and exactly one which contains a — 1; thus, on neighbourhoods of a and a — 1, 
the invariant density has the simple form p^j^x) = x ^ B , 

As an example of such numerical evidence we report a numerical simulation 
of the invariant density for some values of a in the interval [ ^|~ 3 , v/ ^~ 1 ] where 
a matching of type (2, 3) occurs. We fit the invariant density with the function 
A + /(x + B + ) on the interval [max{S , },a] and with the function A-/(x + B_) 
on [a — 1, minis'}]. 



2G 





a = 0.310 


a = 0.320 


a = ^ 


a = 0.338 


a = 0.350 


a = 0.360 


A+ 


1.76114 
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1.78963 


1.81981 


1.84658 


B+ 


1.64768 
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1.62374 


1.62987 
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1.65138 


A- 


1.77289 


1.78874 
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1.84562 


1.85959 


B_ 


2.66097 


2.66081 


2.66583 


2.66751 


2.66915 


2.6658 



Moreover, from these numerical data it is apparent that the leftmost branch 
of hyperbola is nothing else that a translation by 1 of the rightmost one (i.e. 
A + = A_, B_ = B + + 1). 
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Figure 12: Invariant density for a = 0.338 



5.4 Comparison with the entropy 

If / C [0, 1] is a matching interval, the knowledge of the invariant density for one 
single value of a E I plus eq. Q allows us to recover the entropy in the whole 
interval. Let a belong to an interval where a matching of type (k±, ^2) occurs 
and suppose, according to the previous conjecture, that on [as, a] the invariant 
density has the form 

for some A,B 6l and x = max{T™(a), 1 < n < ki}U{T™(a-l), l<m< k 2 }. 
Then by Q, for x < a sufficiently close to a 

m = h[a) (7) 

We think that the entropy has in general such form for values of a where a 
matching occurs. 
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Let us consider the particular case of the interval [0.295,0.3042]. In the 
region to the right of the big central plateau (i.e. for a > ~ 3 + v/ ^ ) the behaviour 
of entropy looks approximately linearly increasing, as conjectured in [9], sect. 3. 
We will provide numerical evidence it actually has the logarithmic form given by 
equation ([7| on the interval [ v ^~ 3 , ;p -]■ To test this hypothesis, we proceed 
as follows: 

1. We fit the data of the invariant density for a = 0.338, obtaining the 
constants A + and B + which refer to the rightmost branch of hyperbola 
(the data are already in the previous table). 

2. We fit the data of the entropy already calculated (relative to the window 
[0.30277,0.3042]) with the function (0. We assume A+ and B+ as given 
constants and we look for the best possible value of h(a) (which we did not 
have from previous computations). The result given is h{a) = 3.28311. 
In the figure we plot the obtained function in the known window, as well 
as a linear fit. In this interval, the difference between the two functions is 



negligible. (Figure 13 1 



In order to really distinguish between linear and logarithmic behaviour 
of the entropy, we computed some more numerical data for the entropy 
far away to the right but in the same matching interval. In this region 
the linear and logarithmic plots are clearly distinguishable, and the new 
points seem to perfectly agree with the logarithmic formulap^j (Figure 14 1 




Figure 13: linear vs logarithmic fit, Figure 14: linear vs logarithmic fit, 0.35 < 
0.3028 < a < 0.3042 a < 0.366 

Notice these data agree with eq. [7] also for x > a, which is equivalent to 
sa y Pa{x) — B ^ +x for a; in a right neighbourhood of a — 1. 

Let us remark that the new values computed are just a few, but are more accurate than 
those in the interval [0.30277, 0.3042] since we used the package CLN a C++ library to perform 
computations in arbitrary precision 
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6 Appendix 

In this appendix we give the proof of two simple results which are of some 
relevance for the issues discussed in this paper. 

Proposition 6.1. If xq is a quadratic surd then Xq is a preperiodic point for 
T a ,ae[0,l]. 

For a = 1 this is the well known Lagrange Theorem, and this statement is 
known to be true for a = and a G [1/2, 1] j8]. Since we did not find a reference 
containing a simple proof of this fact for all a G [0, 1] we sketch it here, in 
few lines: this proof follows closely the classical proof of Lagrange Theorem for 
regular continued fractions given by [5] which relies on approximation properties 
of con ver gents, therefore it works for a > 0. 

If Xq is a quadratic surd then Fo^o) = for some Fq(x) := Aqx 2 + Bqx + 
Co quadratic polynomial with integer coefficients. On the other hand, sinccp^] 
x o = P q lZ\ X x lX P q l ' setting F n (x) := F ( ^-^+^ )(g ra _ 1 x + q n ) 2 , we get that 
F n (x n ) = Fo(xo) = 0. 

Moreover F n (x) = A n x 2 + B n x + C n with 

A n = F Q {p n _ 1 /q n - 1 )ql_ 1 , C n = F ( Pn / 'q n )ql, B 2 — AA n C n = Bq — AA C . 

(8) 

Both A n , B n are bounded since: \F a (p n /q n )\ = \F a (p n /q n ) - F (x )\ = 
l^olOllf 2 " — x o\ < j^2-; moreover from the last equation in ^ it follows that 
B n are bounded as well. 

Proposition 6.2. The variance cr 2 (a) is constant for a G — 1, (v5 — l)/2]. 

This result relies on the fact that for all a G [\/2 — l, (v5 — 1)/2] the maps T a 
have natural extensions T a which are all isomorphic to T\/2- In the following 
we shall prove the claim for a £ [y2 — 1,1/2] and we shall write T\ instead 
of T a and T2 instead of T\/ 2 - So Tj : Ij — > 7y, (j = 1,2) are 1-dimensional 
map with invariant measure [mj; Tj : Ij — > Ij, (j = 1,2) are the corresponding 
2-dimcnsional representations of the natural extension with invariant measure 
fij, and <fr : I\ —> I2 is the (measurable) isomorphism 

$ o fi = f 2 o = /i 2 

First let us point out (see [12] pg 1222-1223) that <!> is almost everywhere dif- 
ferentiate and has a diagonal differential; moreover Tj are almost everywhere 
differcntiablc as well and have triangular differential. Therefore 

d^Txix^dT^^^) = dT2\q>{ x ,y) d^(x. v ) (9) 

and it ie easy to check that, setting IJ the first component of Tj, a scalar 
analogue holds as well 



d<I> x dTf _ dT% d$ x 



(10) 



To simplify notations we shall write p n , q n instead of p n ,a, <?n,a- 
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So we get that, for all k, 



log 



off 
dx 



log 



dx 



o <f> 



log 



d$ x 
dx 



log 



da; 



Since Tf is /2i-measure preserving J j log 
taking into account that = //2 we get 



log 



dx 



oTi 



oTi 



d/ii = 0; so, 





<9Tf 










/ log 


9a; 




/ log 

Jii 


ox 


djli 



log 



Let us define g\ := log 



dT' 



dx 
N-l 



and 52 : = log 



dTZ 



dp.2 



dx 



dx 

(so that g\djli 



Ji 2 92djXi = 0) and S^g := X)fc=o 9°^! we easily see that 

d$ x 



S^fli = 5^sio$log 



<9a; 



log 



which means that Sjj 1 gi and Sjf g2 ° differ by a coboundary. 
Lemma 6.3. Lei u, u 6e iwo observables such that 

1. lim^+oo J(^§) 2 dfi = I e R; 

2. u — v + (f — foT) for some f E L 2 . 
Then 

lim / ( N r — ) 2 dri = lim / (—^=) 2 du. 



The lemma implies 



lim 

iV— >+oo 



/l 



dill = lim 



/2 



N 



(11) 



9Tf 



This information can be translated back to the original systems: since -gj- \ ^ x , y ) 
T l( x )> TErl(*,i/) = ^C 21 ) if we define 

Gi :=log|7t(a:)|- / log 17^)1^ 



G 2 = log|T^(x) 



/2 



log|^(a;)|dM 2 



we get gi(x,y) — G\(x) and g2{x,y) = G 2 (a;); therefore Sj^gi — Sj^Gi and 



^jv 52 = S^G2- Finally, by equation (111}, we get 



lim 

AT— »+oo 



Sj}Gi 



N 



dii\ = lim 

JV->+oo 



JV ^2 



^2 
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6.1 Tables 



fei k 2 ) 
3 9) 

2 8 

3 8 

2 7 

3 8 
3 7 
47 

2 6 

3 8 
3 7 

3 6 

4 6 
9 7 

5 6 
9 7 
10 

2 5 

3 8 
3 7 

3 6 

4 6 
9 7 
3 5 
1 6 
7 6 
1 5 



7) 



7.69e-4 
3.68c-3 
l.llc-3 
5.44c-3 
6.98c-4 
1.69c-3 
8.12c-4 
8.54e-3 
6.06c-4 
1.12c-3 
2.76c-3 
1.34c-3 
2.38e-5 
7.91c-4 
2.25c-5 
1.54c-5 
1.45c-2 
6.57c-4 
1.06c-3 
1.98c-3 
7.42c-4 
1.03c-5 
4.94c-3 
8.44c-4 
l.llc-4 
2.45e-3 



-8+^82 



9 

-4 + s/Vt 

-7+v / 65 



-7+^53 
2 



-6+5^2 
7 

-17+v / 445 



26 

-3+ v/10 

-ll + y/145 
6 

-8 + ^82 



6 

5 + ^37 



-7+y / 82 
11 ' 

-51+13 V29 



100 
-3 + ^145 



53 ■ V5185 
99 ' 

-127+V30629 



250 
-5 + ^29 



10 

-9+yloT 



-13+V229 
10 



14 ' 
-81 + V13226 



155 
-4+y/26 



5 




-10+V 


145 


9 




-25+v 


1297 


48 




-11+V 


JJ 7 ) 



a+) 

-2+y^ 

-7+V77 \ 
14 ) 

-7+3 Vf\ 



6 



-9+2 \/30 
13 



(fcl fc 2 ) 



-3+2 V3 N 
3 



-5+3 Vs A 

10 ; 

-10+2 y/42 N 
-15 + ^357 ^ 



22 

-5+^35 ^ 

-15 + y / 285" '\ 
10 ) 

-117+V15621 N 
42 

-10+2 V30 \ 
5 J 

-30+4 Vtie A 
13 



26 ) 

-io+vISs A 
19 ) 

-4+y3o ^ 

^) 

-7+V69~ \ 

-187+3 y/iolig ^ 
82 , 

-19+3 y/gfl A 
26 j 



-29+V1023 A 

13 ; 

-3+2 V3~\ 



(8 6) 


6.42c-5 


(5 5) 


1.46c-3 


(2 4) 


2.77c-2 


(3 6) 


2.1c-3 


(4 6) 


7.02c-4 


(3 5) 


3.97c-3 


(4 6) 


5.77c-4 


(4 5) 


1.51c-3 


(5 5) 


7.88c-4 


(3 4) 


1.02c-2 


(4 6) 


8.86c-4 


(4 5) 


1.78c-3 


(5 5) 


7.09c-4 


(8 6) 


2.73c-5 


(4 4) 


5.24c-3 


(8 6) 


2.73c-5 


(2 3) 


6.32c-2 


(4 6) 


6.9e-4 


(4 5) 


1.72c-3 


(3 4) 


9.87c-3 


(4 5) 


1.45c-3 


(4 4) 


3.82c-3 


(5 5) 


6.75c-4 


(3 3) 


2.68c-2 


(2 2) 


2.04c-l 


(2 1) 


3.82e-l 



64 

-7+ylOl 
13 

f-2 + V5 

-7+y/65 
4 



' -5+^37 



4 




-13+v 


r 2S7 


11 




-15+v 




22 




-10+v 


22B 


18 




-3+vT7 


4 




-11+v 


170 


7 




-15+v 




14 




-11 + v 


/ 257 


17 




-54+v 


=7^7 



-i+VTf 

7 

-54+V70 
101 

-3+VT3 



13 

-2+vll) 
3 ' 

-l + \/2 

-l+v^ 



34 

-2 + Vtj 

-3+v^ A 

6 J 

-6+4 y/g A 
11 ) 



-9+V165 \ 
14 ) 

-2+2 s/2 \ 
3 ) 



-8+3 yTI \ 
7 ) 
-23+5 y / 29~ '\ 

14 ; 

-3+V15 A 
-19+3^ ^ 

34 ; 

-7+3 yTl A 
10 J 

-6+2 

5 ) 



-■i+s/Vi \ 
2 J 
-127+7 v/453^ 



-1 + v^ 



-13+\/290 


-23 + v / 1365 


11 1 


38 


-15+v / 533 
22 


-4+V37A 
4 J 


-7+V/85 
6 


-3+2 


-9+VT45 


-8+2 v'SA 


8 ' 


13 ; 


-5 + \/65 


-ll + v^22l\ 


8 > 


10 J 


-13+5 ST3 


-2+^^ 



-1+V2) 



A sample of matching intervals found as in section [4] 
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(fci £2) size 
(257 257) 5.43e-201 
(129 129) 7.27e-101 



(a- 



(65 65) 

(33 33) 

(17 17) 

(9 9) 
(5 5) 
(3 3) 
(2 2) 



7.98e-51 

8.81c-26 

2.78c-13 
5.2c-7 
6.75e-4 
2.68c-2 
2.04e-l 



{ -1051803916417 
+ 5 \Al0424870216034832616745 }/ 
1576491320449 



-1 + 



s/3152!)h2n 10!) 
128045 



-433+\/4i<7.v~,7 
649 

-13+5^13 
13 ' 

-2+vTo 
3 ' 

-1 + V2 , 



, V31529826409 
' 128045 

-433+V467857 
649 

-13+5 y/13 



13 



-2+VT0 



V2) 



-1+V5 
2 



A chain of adjacent matching intervals (see section 4.2 1 
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